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The model of community isolation was extended to the case when indi- 
viduals are randomly placed at nodes of hierarchical modular networks. It 
was shown that the average number of blocked nodes (individuals) increases 
in time as a power function, with the exponent depending on network pa- 
rameters. The distribution of time when the first isolated cluster appears 
is unimodal, non-gaussian. The developed analytical approach is in a good 
agreement with the simulation data. 

PACS numbers: 89.75.Hc, 02.50.-r, 89.75.-k, 89.75.Da, 89.75.Fb 

1. INTRODUCTION 

Recently, hierarchical systems have been attracting attention of scien- 
tists working on complex networks [U O [3l Ul |5] . In fact many real networks 
are hierarchically organized, e.g. WWW network, actor network, or the 
semantic web pQ. Dynamics at such networks can be qualitatively and 
quantitatively different from that at regular lattices (see O O H] ) . 

The Ising model at a network with a hierarchical topology was studied by 
Komosa and Holyst [2j. The analyzed parameters were, among others, mag- 
netization, magnetic susceptibility, critical temperature and correlations of 
magnetization between different hierarchies. It was shown that the critical 
temperature is a power function of the network size and of the ratio ^j^, 
where k stands for a node degree. 

Opinion formation in hierarchical organizations was studied by Laguna 
et al. [3]. Agents, belonging to various authority strata, try to influence 
others opinions. The probability that an opinion of an agent of a certain 
authority prevails in the community depends on the size distribution of 
the authority strata. Phase diagrams can be obtained, where each phase 
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corresponds to a distinct dominant stratum (or a sequence of the strata, 
with the decreasing probabihty of prevaihng). 

Fashion phenomena at hierarchical networks were studied by Galam and 
Vignes Interactions were imposed between social groups at different 
levels of hierarchy. A renormalization group approach was used to find 
the optimal investment level of the producer and to assess the influence of 
counterfeits on the probability of a new product success. 

One of fundamental topics in social dynamics are conflict situations and 
many different sociophysics approaches [6l[71[8] or prisoner's dilemma-type 
games [9] have been proposed. Recently a simple model of communities 
isolation has been introduced by Sienkiewicz and Holyst [10] . The model 
can describe such various issues as strategy at battlefields or formation of 
cultures. The idea behind this model is similar to the game of Go and it 
takes into account a natural leaning of people to avoid being surrounded by 
members of another (potentially hostile) community 

In this paper we extend the model of communities isolation studied for 
chains, hypercubic, random and scale-free networks |10[ [T2] to hierarchical 
networks proposed by [1] . 

2. HIERARCHICAL NETWORKS 

The model of hierarchical networks was proposed by Ravasz and Barabasi 
[1] and modified by Suchecki and Holyst [5j. Such networks possess 3 pa- 
rameters determining their structure: 

• The degree of hierarchy /i E N U {0} 

• The distribution PM{m), where m G N, determining number of nodes 
at each level of hierarchy (in particular, the size of the cliques at the 
lowest level of hierarchy is m + \) 

• The parameter determining the density of edges p G [0,1] 

Two models (referred to as PI and PD models) were analyzed, which differ 
in the density of edges. Each network has a central node, referred to as a 
center of hierarchy. A network of hierarchy /i = is a complete graph of 
size m + 1 (m is a random number, chosen with probability Pa/ (m)). The 
center of hierarchy, due to the symmetry, is an arbitrary node. In order to 
construct a network of hierarchy h > 0, one has to construct m + 1 subnet- 
works of hierarchy h — 1 and choose one of them — its center of hierarchy 
becomes a center of hierarchy v of the whole network. Afterwards, new con- 
nections (edges) are created: for each node w of m remaining subnetworks 
a connection (edge) {v, w) is created with probability p (in case of the PI 
model) or p'^ (in case of the PD model). Sample networks created this way 



3 




are presented in fig. [TJ Let us stress that the subnetworks do not have to 
be connected, especiaUy if p is smah. 

Some basic properties of such networks can be concluded from the con- 
struction algorithm: 

• For h E {0, 1}, as well as for p G {0, 1}, models PI i PD are equivalent. 

• For p = the network consists of isolated cliques of size m + 1 [m — 
random variable). 

• For PM{m') = Sm',m the number of nodes (vertices) of the network 
equals N = {m + 1)'^+-'^. 

Periodic oscillations in degree distribution of such networks can be observed 
in log-log scale. The period, the amplitude and the shape of the peaks 
depend on the parameters of the network [5J. 

In this paper only the case with PM{m') = 5m', m {rn = const) was 
considered, which corresponds to the original Ravasz and Barabasi model 

pa- 

3. BASIC ISOLATION MODEL 

The model of communities isolation was proposed by Sienkiewicz and 
Holyst [TD]. The rules are similar to those of the game of Go. A number 
of communities compete with each other, settling nodes of a network. In 
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each step a random empty node is chosen. It is then settled by a member of 
randomly chosen community. A cluster of nodes occupied by one commu- 
nity becomes blocked when it gets surrounded by another community. The 
surrounded nodes are no more active in the game, i.e. they can not take 
part in surrounding other communities. 

The case of communities competing at a chain was analyzed in [lOj . 
Two functions describing the evolution were studied: the average number 
of blocked nodes over time and a mean critical time, i.e. the moment, when 
the first blocked cluster appears. In [T^] the influence of external bias was 
considered when settling rates of competing communities are different. 

In this paper the case of two competing communities at PI and PD hier- 
archical networks is considered. Two parameters are analyzed: the average 
number of blocked nodes Z{t) and the critical time distribution Pr{tc)- 



4. NUMBER OF BLOCKED NODES OVER TIME 

4.1. Case p = 

For p = models PI are PD equivalent. The network consists of ^ 



m+l 

isolated cliques of size m + l. In such case 

Z(0) = , 

z(t + i) = z(t) + E.=iip., ^ ' 

where pi — probability that in the {t + l)th step i nodes will be blocked, 

t \ I m 



After short algebra we obtain 

^(0) = , 
Z(t + l) = Z(t) + f (^)™. 

The solution of this recursive equation is a (m -|- l)th degree polynomial, 
which can be approximated by substituting the sum with the integral: 

™ +m+l 
= _ 

m+l 2iV'" ^ ^ 

As one can see, Z{t) is a power function. The exponent (3 depends only 
on the m parameter, f3 = m + 1. 
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Case p = l 

In this case models PI and PD are also equivalent. For networks of 
hierarchy h = 1: 
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^mpopT{l + {m + l)pi), (5) 



where pi is a reduced density: 

for t > i 



, . , for t < i 

Af-i 



For networks of higher hierarchies, /i > 1, a recursive equation well ap- 
proximating Z^^\t) can be derived. The idea behind the formulas is that 
a clique can only be blocked if all the nodes of higher hierarchies neighbor- 
ing with it are filled. Therefore Z^^\t) = if the center of hierarchy of 
the network (which neighbors with all the other nodes) is empty. In the 
opposite case, Z^'^^t) depends on Z^^^^\t), which describes each of m -|- 1 
subnetworks. 



Ziy) = ^mp,p-l,il + im + l)p,+i) 



Zt\t) = zt'\t) + \mp,Z\l-'\t) 

+ lp^P^+,P^::'^"'\{m+l)'^ + l) 

Z^^\t) = Z^^\t) 

This equation can only be solved numerically. The solutions are presented 
in fig. [2] and [3l It can be noticed, that within a wide range of time t, Z^'^) 
can be with reasonable accuracy approximated with a power function 

zW(t)oct'^, (8) 

where the exponents /3 are higher than in the case of p = and they are 
close to m + A. 



4-3. General case 

An analytical approximation of Z^^\t) for networks with higher hier- 
archies {h > 1) when the parameter p is different from zero and one is far 
more difficult. Instead of searching for such a formula, an alternative ap- 
proach was chosen. It was assumed that Z^^\t) can be estimated from the 
proportion 

log Zl%{t)- log Z(^Ht) ^ f(p,h) 
log ZW{t) - log 4=1 (*) ^ 1 - /(^'^ ^) ' 
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where f{p, h) € [0, 1] should be an increasing function of p which, while 
not being too complicated, would give a reasonable approximation for the 
widest possible ranges of p and h. It turned out that in the case of the 
PI model, choosing f{p,h) = p results in a good agreement of the Z^^\t) 

function with simulation data. For the PD model, /(p, h) = pa is a good 
choice. 



5. CRITICAL TIME DISTRIBUTION 

5.1. Case p = 

As it was previously mentioned, in the case of p = the network consists 
of = (m + 1)'* isolated cliques of m + 1 nodes. In order to find the 
distribution of critical time (time, when the first blocked cluster appears), 
one has to consider the probability that at time t there are no blocked nodes 
yet. It means that at time t the only completely filled cliques are those filled 
with members of one community, which leads to the formula 



Pr{tc > t) 



m+l^ 



N 
m+1 



(10) 



where a = 1 — 2 The cumulative critical time distribution can be im- 
mediately obtained 



Pr{te <t) = l- 




as well as the critical time distribution in the approximation of continuous 
time: 



Pr{tc = t) = Pr{tc <t)- Pr{tc <t-l) 
« j^Prit, < t) 

N 1 

(/ J. \ m + l / + \ '"^ 



7 



I h = 3, m = 3, model P1 | 
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Fig. 2. (color online) Average number of blocked nodes, Z{t), for various networks 
of PI model. Symbols correspond to data from computer simulations. Lines show 
analytical approximations (eq. (jH), ([7]) and Left side — linear scale, right 

side — log-log scale. 
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Fig. 3. (color online) Average number of blocked nodes, Z{t), for various networks 
of PD model. Symbols correspond to data from computer simulations. Lines show 
analytical approximations (eq. (jH), ([7]) and Left side — linear scale, right 

side — log-log scale. 
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The mean critical time can be also calculated analytically: 



(tc) = tPr{tc = t)dt^ t—Pr{tc< 



t)dt 



N i_ , „ 

-a "i+i B 



m + 1 



1 

+ 1, 



m + 1 



m + 1 



B 2 



N 



m + 1 



+ 1, 



1 



m + 1 



(13) 



where B{a, h) = t°- ^(1 — if ^ dt (Euler beta function) and B{x; a, h) = 
Jq t°--^{l - tf-^ dt (incomplete Euler beta function). 



5.2. Case p=l 
For networks of hierarchy h = 

^ \ m+1 



Pr{tc >t) = l-a[- 



For networks with hierarchy h = 1 

Pr{U >t) = l-p^ + po(l - pT) (1 - ap^')' 



(14) 



(15) 



For networks with any degree of hierarchy, /i > 0, a recursive formula for 
the cumulative critical time distribution can be expressed as 



r Fi'\t)=apr' 



Ft\t) = p,- p,{l - pT+,) ( n (1 - Fll\{t)) 

\d=0 

[ PA''\t,<t) = F^''\t). 



h-l 



(16) 



The mean critical time can be obtained by numerical integration of 

Pr^''\tc < t): 



(t^) = N - Pr^^\tc<t)dt. 
Jo 



(17) 



6. DISCUSSION AND CONCLUSIONS 

6.1. Number of blocked nodes Z(t) 

In all cases the function Z(t), defined as the average number of blocked 
nodes at time t, can be approximated with a high accuracy by a power 
function 

Z{t) oc t^, (18) 
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Fig. 4. (color online) Critical time distribution for networks of hierarchy h = 2 (left) 
and h = 3 (right), with p = (top) and p = 1 (bottom). Symbols — simulated 
data, smooth lines — analytical approximation (eq. and ([TCI)). 



The /3 exponent depends on the parameters of the network. For d-dimensional 
hypercubic networks (including the 1-dimensional ones, i.e. chains) /3 = 
2d — 1. For modular hierarchical networks /3 depends mainly on m and p 
parameters, i.e. on the sizes of basic cliques at the lowest hierarchy and on 
the density of inter-clique connections. The dependence on the degree of 
hierarchy h (and on the network size) is weak, what can be explained by the 
fact, that increasing the degree of hierarchy /i is a process similar to system 
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Fig. 5. (color online) Mean critical time for various networks. Symbols correspond 
to data from the simulations, lines — to the analytical approximations (eq. (jl3|) 
and (HZl)). 

rescaling, therefore 

Z('^+I)(p) « (m + = (m + 1)C7/ = C'tl^. (19) 

For p = 0, the parameter /3 can be analytically found: /3 = m + 1. 
The result is in agreement with the simulated data. Increasing the den- 
sity of connections (the p parameter) leads to the increase of /3 — up to 
approximately m + 4 for p = 1. 

There is an important distinction in the way Z{t) was approximated 
for hypercubic and hierarchical networks. For hypercubic networks, the 
number of isolated nodes was calculated using the following approxima- 
tion: all blocked nodes were blocked alone, i.e. they do not neighbor with 
other blocked nodes (of the same community). Although this approxima- 
tion might seem coarse, resulting analytical predictions turned out to be 
in quite good agreement with simulated data \10\ [T2] . For modular hierar- 
chical networks, such an approximation would not be reasonable. Because 
of the fact that at the lowest level of hierarchy such networks consist of 
cliques of m + 1 nodes, the most probable are situations when ^^^^^ nodes 
are simultaneously blocked. 

6.2. Critical time tc 

The second analyzed parameter was critical time tc, i.e. the moment, 
when the first isolated cluster appears. It is a random variable. The critical 
time distribution Pr{tc) was studied, as well as mean critical time {tc)- More 
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precisely, a critical density (or a critical relative time) 

Pc^'^ (20) 

was often shown so networks with different parameters could be easily com- 
pared. 

The Pr{pc) distribution is always unimodal. The mode (arg max Pr(pc)) 
decreases with the increase of h and the standard deviation cr{pc) decreases 
with m. 

For p = it was possible to find the analytical formula for both Pr(tc) 
and (tc). The distribution Pr{tc) is a polynomial of degree m{m + l)((m + 
1)'* — 1) (see eq. \T2\i and the average (tc) is a scaled difference of two Euler 
beta functions (see eq. [T3]l . The average {pc) decreases with h and for a 
fixed h it reaches a minimum for m ~ 2 (see fig. [5]). 

For p = 1 the distribution Pr{pc) reaches a constant, non-zero value for 
/?c G [1 — e, 1] (for h < 3, e 0.1), which means that processes when blocked 
clusters firstly appear at the very end of the evolution are not unlikely. 

The values of (pc) can be compared with those obtained for hypercubic 
networks. Similar trends can be observed in hypercubic and hierarchical 
networks: (pc) decreases with the network size N and increases with the 
average degree. However, for modular hierarchical networks the dependence 
of (pc) on the average degree (which equals m for p = and rises with p) 
is very weak in comparison to hypercubic networks. Typical values of {pc) 
for hierarchical networks correspond to the ones obtained for two- or three- 
dimensional networks, even for m ^ 3. 

Acknowledgments 

The authors acknowledge support from the European COST Action 
MP0801 Physics of Competition and Conflicts and from the Polish Min- 
istry of Science Grant No. 578/N-COST/2009/0. 



REFERENCES 

[1] E. Ravasz, and A.-L. Barabasi, Phys. Rev. E, 67, 026112 (2003). 
[2] S. Komosa and J. A. Holyst, Ising model at hierarchical network, to be pub- 
lished. 

[3] M.F. Laguna, S. Risau Gusman, G. Abramson, S. Gonalves and J.R. Iglesias, 
Physica A 351, 580 (2005). 

[4] S. Galam and A. Vignes, Physica A 351, 605 (2005). 

[5] K. Suchecki and J. A. Holyst, Acta Phys. Pol. B 36, 2499 (2005). 



13 



[6] I. Dornic, H. Chat, J. Chave, and H. Hinrichsen, Phys. Rev. Lett. 87, 045701 
(2001). 

[7] G. Deffuant, F. Amblard, G. Weisbuch, and T. Faure JASSS 5 (2002). 
[8] S. Galam, Eur. Phys. J. B 25, 403 (2002). 

[9] S. Lee, P. Holme and Z.-X. Wu Phys. Rev. Lett 106, 028702 (2011). 
[10] J. Sienkiewicz and J. A. Hoiyst, Phys. Rev. E 80, 036103 (2009). 
[11] T.C. Schelling, J. of Math. Soc. 1, 143 (1971). 

[12] J. Sienkiewicz, G. Siudem, J. A. Hoiyst, Phys. Rev. E 82, 057101 (2010). 



